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^ Abstract. A family of algebraic surfaces with many nondegenerate real singularities is 

^ introduced with the help of a construction, which has been used in previous works for the 

' ^ ^ generation of substitution tilings. 
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1. Introduction 



In 1819 a quartic wave surface was discovered by Fresnel in the study of crystal optics. Later 
Kummer observed that the surface had the maximum possible number of nodes of a quartic in 
the complex projective space P^(C) and he constructed a surface with only real nodes. Since 

^^ then the subject of algebraic surfaces has a rich history with plenty of achievements. 

"^ A construction of surfaces with degree m in P^(C) having only nondegenerate singular points, 

; or double points, was introduced by Chmutov in [4J. It was based on his previous idea of using 

^ Chebyshev polynomials and on the Givental methods for generating cubics [1]. 

^_^ In [3] it has been shown that the Chmutov construction can be adapted to give surfaces with 

1— I only real singularities. In C^ a node or Ai singularity can be written in the form Wi+W2+w'^ = 

^ in some local coordinates. In R^ there are two types of nodes, conical nodes (x^ + y^ — z^ — 0) 

and solitary nodes (x^ + ^^ + z^ = 0) and the construction produces conical nodes. The main 
idea consists in defining real folding polynomials F^^^{x^y) G R[x,y] associated to the root 
C^ lattice A2 (which are generalizations of Chebyshev polynomials, related with the root lattice 

Ai). In the cases with degrees 6,7,8,10,12, there are higher lower bounds for the maximum 
number of real nodes /xr on a surface of degree m. A 65-nodal sextic and a 345-nodal decic 
invariant under the symmetry group of the icosahedron were constructed by Barth [2j. Other 
surfaces with the higher lower bounds known are Labs's 99-nodal septic, Endrass's 168-nodal 
octic, and the surface of degree 12 with /iR = 600 nodes, invariant under the reflection group of 
the regular four-dimensional 600-cell, introduced by Sarti O [HI [121 [13] . The results in ^ show 
that all known lower bounds for the maximum number of nodes on a surface of degree 772, can 
be attained with only real singularities. 

In this article we consider a construction motivated by recent work on tilings generated by 
substitution or inflation rules. In particular, two subfamilies of surfaces with degree divisible by 
3 are obtained. One is related to the dihedral group D-^ and the other to the cyclic group C3. 
We use Mathematica [15] and Surfer [6l [13] computing and geometric visualization tools. 




Figure 1. Simple arrangements T^d — {D)odd with seven, nine, eleven and twelve lines as 
subarrangements of the systems of lines (D). li m — 'iq then Yj^ has i?3 symmetry. 



2. The construction of surfaces with degree m 

Arrangements of lines are labeled sets of lines not all passing through the same point. A vertex, 
which is the intersection of two or more lines, is called ordinary if only two lines meet. If all the 
vertices are ordinary, then the arrangement is said to be simple [lOj. 

By a d-system we mean a simplicial arrangement of straight lines in d orientations. There are 
four types of geometric constructions, denoted (A), (B), (C) and (D) in [^ which are the basis 
for the generation of tilings with arbitrarily high symmetry. They are of interest in the study of 
quasicrystals and other fields [8^^ ^9] . For our purposes we may consider just two types because 
the systems with d odd are included in the systems with d even: (B) is included as a subsystem 
of (A), and (D) is included in (C). In this work we denote the two basic constructions by (D) and 
(C) which correspond to (A) and (C) in |8j. We donote Sy = sin{vTi /d)^ty = tan{vTi/d). The 
following sets of lines in the xy plane define, for each d, up to mirror images, one or two finite 
patterns made with triangles having edge lengths Sy^ and containing the necessary information 
for the derivation of substitution or inflation rules that generate aperiodic planar tilings: 



where L'^^{x^ y) := y — xt^ — F^/ n and the term F^/ x, X=C,D is defined as follows 



(1) 



• (D) d = 2m, (m = 3,4,...) 

For V — 1,2,.. .,771 — 1, the index a{v) is defined as a{v) — v {v 

a{v) — m — v {u — \m/2\ + 1, ...,7n — 1) and F^/ n is defined by F^/ n 

For z/ = 771+ 1,771 + 2, ..., d — 1, we have a{v) — v — m (z/ = 77i + l, ..., 77i+ [77i/2j), a{'u) — d 

{y^m^ [m/2j + 1, ..., d-\) and F^^^ '.-- 
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771+1, ..., 2771 — 1 we 

m+ l,...,m+ [(m- l)/2j) 



and a{v) = 2m-l-z/(z/ = m+ [(m - 1)/2J + 1, ..., 2m - 1). 



The 2771-systems (X) are the union of two subsystems: (X) 



even Contains L^^2m{^^y) 







for z/ even and {X)odd has L^2m(^5^) = for z/ odd. We consider, in what follows, a 
kind of constructive duality between substitution tilings and surfaces with many singularities. 
The subsystems {D)eveni {C)odd are on the basis of the general geometric constructions for 
substitution tilings with 771-fold symmetry (see [9J for the subsystems corresponding to 771 = 



5, 10), which have been studied in [8] and previous works [7J. The subsystems T^n — {D)oddi ^C — 
(C)even ^ire simple arrangements that wih be used as the starting point of the following 
construction of surfaces with degree m having many real singularities. We define two families 
of polynomials associated with each simple arrangement: 
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Figure 2. Simple arrangements T^c — {C)even with nine, twelve, fifteen and eighteen lines as 
subarrangements of the systems of lines (C). They have C3 symmetry. 



Let Tjn(z) be the Chebyshev polynomials of the first kind generated by the recurrence relations 
Tm^i{z)-2zTm{z)+Tm-i{z) = 0, m > withro(z) = l,Ti(z) = z (Fig.3). The real projective 
surfaces of degree m associated with S = Ti]j^ S^ are defined by 



mTi ^ 



Pms{x,y,z) := Jms{x,y) + -—-{Tm{z) + 1) G R[x,y,z] 



(3) 



where —Xm^ is the unique critical value (minimum) of JmJ:{x, y) in the interior of the triangular 
regions of the simple arrangement S. In Fig. 4, which corresponds to ttt. = 9, we have shown in 
black the regions where JmJ:{x,y) is negative and the minima are located inside the bounded 
triangular zones. The three different triangular shapes in T^^ and the seven triangular shapes 
in T^c form two sets of prototiles for quasicrystal substitution tilings as shown in [7J. The 
polynomials JmJ:(x^y) have critical points with only three critical values 0, —XmJ: and A^j], 
therefore the surface Pjnj:(x^y^z) = is singular only at the points where the sum of the 
critical values of Jm^{x,y) and ^^(T^(z) + 1) is zero: either both are 0, or the first is 



—XmJ: and the second A^^. The arrangements S^ are equivalent to the simple arrangements 
corresponding to the real folding polynomials F-£i^{x^y) in [3J. In fact it is possible to 
relate JmY^oi^^y) ^^^ ^Rm(^'^) ^^ introducing a change of variables and a scaling factor: 



if m is even 



(m+2)7r . (m — 2)7r 
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^^^^^^^ ^ if 771 is odd. In Fig.l we have shown the cases tti = 7, 9, 11, 12 
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and 6^ - a^ - 3 • v_^.^3. 
for Ej). When m = 3g, the symmetry group of E^ is the dihedral group D^,. 
The surfaces PwEd (a^i J/i -2^) = ^^d |3j 



,^2 



^^2 



1 



Chm^^{x,y,z) := F^^(x,y) + -(r^(z) + 1) = 



(4) 



have the same number of nodes. The number of points of JmT.{x^y) with critical value zero is 
equal to the number of vertices in the simple arrangements, namely (^) . The number of triangles 
in Y^D is \m? — m ii m — mod 3 and ^m? — tti + | otherwise. On the other hand Tjn{z) have 
[yj points with critical value —1 and [^^^^J points with critical value 1, hence the number of 
real nodes of PmT.D (^^ Vi^) ^^ and Chm^^{x, y^z) = is /xr = (^) [ yj + (|m^ - m) [^^J if 
m G Z3 and fin = {'^) [f J + (^m^ - m + |) [^J otherwise. 





Figure 3. The Chebyshev polynomials Tg(z) and Tis{z) 




Figure 4. The polynomials J9^^(x,y) (left) and J9^^(x,y) (right) have negative values within 
the regions in black, and minima inside the bounded black triangles. 



3. The surfaces with degree m — ?>q associated with Sc 



Within Y^c there is one more triangle than in S/^ while the number of vertices is the same, 
therefore the corresponding surfaces for m — ?>q have more singularities. As we have done in 
the previous section, in relation with S^, we can also make a change of variables in Jm^ci^^ v)^ 
which facilitates our analysis. We translate the origin to the barycenter of the unique equilateral 



triangle in S^, having coordinates (c3g, |) with C3g 
corresponding to the minimum in eq.(3) is then 



\/3 I gm((g-l)7r/6g) 



+ 



2sin{7T/6q) 



The parameter 



AsgEc = Jsq^ci^^q^ 



(5) 
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Figure 5. Local miniina with value -1 obtained when Sc for tti = 9 is rotated by /c7r/54. From 
left to right: 7r/54 (prototiles /,^, e, rf),27r/54(^, a), 87r/54(^, c), 147r/54(6,^). 



The symmetry group of the simple arrangements S^ is the cyclic group Cs (Fig. 2), hence for 

+ 1 different 



the study of the critical points of JruT^ci^^v) '^^ enough to consider the ^ g~ ^ 
triangular shapes which are, on the other hand, the prototiles of certain classes of substitution 
tilings ^ i8j . The action of C3 gives an orbit of length 3 except for the central triangle therefore 
we have \im? — tti + 1 minima. A first set of minima (which correspond to the triangles above 
the horizontal line and with an edge on it in Fig. 2) can be obtained by rotating Yjq an angle of 
O^ ~ 6m ^^^ ^^^ remaining ones with rotations by integer multiples of O^. We also change 
the scale in order to have all the local minima equal to —1 and all the local maxima equal to 8. 
We define for m — 'iq 



1 



with 



P^ (x, y, z) := J^(x, y) + -{Tm{z) + 1) G R[x, y, z] 



u = ajnXcosQjn + amysinQm -\-Cm, V = -amXsinQm + amycosQm + z 



(6) 



(7) 



where the coefficients a^ have been given in Section 2. 

The surface of degree 9 corresponding to T^c in eq(3), contains the polynomial 

J9Sc(x,y) = {y - xt2 + 57 + S5){y - xt4 + 57 + 55 + 53 + si){y - xt^ + 57 + s^){y - xts) 

{y - xtio - 59 - S7){y - xti2 - sq - S7 - s^ - ss){y - xtu - 59 - 57 - 55 - 53) 
{y - xtie - 59 - S7)y 

(8) 
and Xqj:^ — ^ 'p^ 13.28. By using the notation in [7J we have three scalene triangular 
prototiles a, 6, c in the 9-system (D)^^^^ and in {C)odd we have four additional triangles: 
one is equilateral {g) and three are isosceles (d^e^f). In H^ (Fig. 4) the triangles above 
the horizontal axis and with an edge on it are, from left to right /, ^,e, d. The tile b has 
a common vertex with / and c shares a vertex with g and another with e. In this case 
J^(x,0) = -1 + 27x2 - 9x^ - 54x^ + 36x^ + 21x^ - 27x'^ + 9x^ - x^ has four minima 
corresponding to f^g^e^d. The minima inside a^c^b are obtained by restriction to the x- 
axis of the polynomials obtained from eq.(7) when O^ is replaced by the rotation angles 
27r/54, 87r/54, 147r/54, respectively (the minimum in the negative x-axis for 27r/54 in Fig. 5 is 








Figure 6. Local maxima with value 8 obtained when S^, for m 
/c = 0, 6, 12 (from left to right). 



9, is rotated by A:7r/54 with 



not -1 and does not correspond to a local extremum inside a prototile). The whole set of 
minima can be obtained with the rotation angles 7r/54 + /tt/S, (2 + 6A:)7r/54 + In/ 3; A;, / = 0, 1, 2 
(Fig. 5). The maxima are located in nine bounded non-triangular regions with three different 
shapes. They can also be reached by rotations of A:7r/54 (see Fig. 6 for the cases with / = 0): 
6fc7r/54 + /tt/S; A:,/ = 0, 1,2. There are (2) = 36 vertices where jQj2^{x^y) has critical value 
and the minima are located inside each of the 19 black triangles in Fig. 4. Now Tg{z) has critical 
value -1 in four points and 1 in also four points, therefore the number of nodes of Pq{x^ ?/, z) = 
is /XR = 220. 

The polynomial JisecI^^I/) ^^^ ^^^ same minimum value — Ai5^^ = — ^ ^ —2.4x10 inside 
the 61 triangles of ^c (Fig.2). Now J^^{x, 0) = -1 + TSx^ - 25x^ - 450x^ + 300x^ + 895x^ - 
945x'^ - 495x^ + 1045x^ - 297x^° - 285x^^ + 260^^^ _ gg^is ^ 15^14 _ ^15 j^^^g ^^^^^ minima with 
value -1, and the minima of the other 14 prototiles are located by rotations of type |^, /c G Z. 



has /xr = 1162. 
2^ and AisEc ^ 4.8x10 . The arrangement He has now 91 triangular 



The surface Pil{x^ y, z 
For m even A^^^ = 
regions (Fig.2), the polynomial Ji^(x,0) = -l + 108x2-36x^-945x^ + 630x^ + 2919x^-3024x'^ 
3366xS + 5720x^ - 4212x^1 + 2457x^2 _^ 378^13 _ lOSSx^^ + 528^1^ - 135^1^ + 18^^^ - x^^ has 8 
minima, the directions for finding the other minima are given by ^, fc G Z, and the number of 
nodes of Pil{x^ y, z) = is /xr = 2105. 



4. The local extrema of J^{x, y) and the number of real nodes of P^{x, y, z) — 0. 



C( 



It is possible to get an alternative formulation which allows to obtain the positions of the extrema, 
having in mind that they are situated in the directions given by ^. For tti = 3g, g = 1, 2, 3, ... 
we define the lines Lk^r^{x, y) = with 



- . . . k27r . kn . k27r 
Lk,m[x, y) :=y - [cos- x)tan- sin- — 

0771 0771 0771 



(9) 



where A: G Z and Lsm,m(x^y) = is interpreted as the line x = —1. The polynomials obtained 
with Lj^^jyi(x^y) are defined as 



i-(-i)" 



m—l 



Jg{x,y) := 3^^(-l)L— J+i f] L6.+i,„(x,y) G R[x,y] 



(10) 



iy=0 



I ^T^~3 I 

uU ^{hnh+2k} (11) 



The non-integer coefficients a in J^(x^y) and J^(x^y) = J^(x^—y) appear only in the 
monomials ax^y^ with k odd. On the other hand, all the monomials in Fj^^(x, y) have integer 
coefficients and, in addition to the relation given in Section 2, we have another one for the folding 
polynomials: F^l^{x, y) = 6 - J^(x, y) - J^(x, y). 

In order to locate the minima we consider the sets of lines M_i = {L^y^r^ — 0, L6i/+2,m — 
0},z/ = 0, 1, ...,771 — 1. The minima are placed at all the points where exactly three lines of 
M_i meet (Fig. 7). We relabel the hnes as Ik^ A: = 1, 2, ...2m, with /i : Lo,m = and the others 
are assigned as they appear in consecutive clockwise orientations. When {//^i H lk2 H lk?>} is 
non-empty, then the subindices are either all even with fcl + A:2 + A:3 = 2m + 4, or all odd with 
fcl + A;2 + A:3 = 2m + 3, modulo 2mi. The set of minima in the x-axis are (we indicate only two 
of the three intersecting lines) : 

m-3 I 
- 2 

The minimum in the central equilateral triangle corresponds to k — q — 1 and the number of 
minima in the x-axis is 1 + L^^T^J- ^^^ remaining minima, connected with intersections with 
odd subindices, are located, up to 3-fold rotation and whenever the subindices are positive, in 

{^3 n /s, ^3 n /t, ..., /3 n hq-?} U {hq-l n /6g-3, hq-1 H hq-^, ..., hq-1 H lQq_^2q-b)} 
{/t n /g, /r n /ii, ..., I7 n /2g-7} U {/6g-5 H /og-T, ^6g-5 H /6g-9, •••, kq-5 H /6g-(2g-g)} 

The last term of the sequence is given by {Iq-i H /g+i} if q is even, or {lq-2 H /g, /g_2 H /g+2} if 

£_2 g-i 2 

g is odd. Now we have 2 X]|=o (^^ + -'-) i^iii™^ if Q is even and 2 Xlj[.lo (^^ ^ ^) ^^ ^ ^^ ^^^• 
Minima appearing in the intersections with even subindices are located in 

{I2 n /4, h n /e, ^2 n /§, ..., h n /2 J u {h n /e, k n /§, ..., ^4 n /2 Ju . . 

u{/6 n /§, ^6 n /lo, ..., /6 n hq} u ... u 1/2^-2 n hq} ^ ' 

which gives Ylk=o ^ poiii^s. Due to the cyclic symmetry, if a minimum appears in l^i H lk2 then 
there are two more situated in lki+2qn H lk2+2qn H l2m+3-ki-k2+2qn^ n = 1, 2. The total amount 
of critical points with value —1 is then 1 + IBxni — Z 

Also the maxima are located at all the points of intersections of three lines ( Fig. 7 ), belonging 
now to the set Ms — {^6z/+4,m = 0},z/ = 0, 1, ...,7ti — 1 which is a TTi-system of type {C)odd- 
In this case the lines are Ik^k — 1,2,...7ti, with /i : ^4^^ = 0. For those sets {Ikiilk2ilk3} 
intersecting in one point we have that their subindices satisfy fcl + A:2 + A:3 = 1, modulo mi. The 
maxima are found, up to 3-fold rotation, in the following YlX=o ^ poi^^^s 

{h n I2, h n /a, h n h, ..., h n /J u {h n h, h n h, ..., h n / Ju . . 

...u{/^_2n/^_i,/^_2n/Ju{/^_in/J ^ ^ 

The remaining ones are located in l}^ij^qnnlj^2+qn^h-ki-k2+qni ^=1,2, when lki^lk2 is included 
in eq(14). Adding up all the points we get ^^^~ ^ maxima with value 8. 

An alternative way of obtaining the number of minima consists in observing that M_i 
is formed by the union of the TTi-systems (D)even and {C)odd (^ mirror reflection of the one 
appearing in the analysis of the maxima) . There are 1 + ^ g~ ^ points of intersection of three 
lines in (D)even ct^d ^ g~ ^ in {C)odd- By taking into account the total number of minima we 
can compute the number of real nodes of P^(x, y, z) = and the result is 

im\.m. , m(m — 3)..m — l. ,^^. 




m=6 



nn=9 



Figure 7. Local minima and maxima of J^(x^y) are placed at all the points where three 
black lines, belonging to the sets M_i and Ms respectively, meet. The simple arrangements are 
represented in blue. 



— ^J the number of nodes of PjnT.^{x^y^z) = or Chm^^{x^y^z) — for 
3g. Instead of the polynomials T(z), we can use J^{z, 0) and we obtain surfaces with the 



which exceeds in [ 
m 



same number of nodes: 



^C 



C( 



Q^{x,y,z) := J^(x,?/) + 



(_l)m+i 



(J^(z,0)-l + (-ir+i2)GR[x,^,z] 



(16) 



The surfaces Q^{x,y^z) = 0,7n = 6,9,12,15,18, visualized with the Surfer program, which 
shows the images within an invisible sphere, are represented in Fig. 8. They have mirror 
partners Q^{x^y^z) = defined as in eq(16) by replacing J^ with J^ ( Q^^x^y^z) — can 
be seen in Fig. 8). We observe that this is not the case for the surfaces Pm^^{x,y,z) = or 
Chm-^^{x^y^z) = 0, because T^d has one {m ^ 3q) or three {m = 3q) mirror symmetry planes 
(Fig.l). 

The polynomials introduced in this work can also be used for the construction of hypersurfaces 
with many real nodes. For instance, in addtion to the Chmutov hypersurfaces in 4-dimensional 
space defined with the folding polynomials [4| we have, in the cases corresponding to tti = 3g, 
new types of pairs of mirror three-folds with affine equations Jml^o, ^i, ^2, ^3) = Jmi^o^xi) — 
J^{x2,xs) = and Jm{xo,xi,X2,xs) = J^{xo,xi) - J^{x2,xs) = 0. We find 3q{q - 1) more 
nodes in these mirror pairs than in the Chmutov hypersurfaces. Also the polynomials J^{x^y) 
can be of interest in the study of singularities of type Aj^j > 1. 
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